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Abstract. The addition relation for the Riemann theta functions and for its limits, which lead to the 
appearance of exponential functions in soliton type equations is discussed. The presented form of addition 
property resolves itself to the factorization of N — tuple product of the shifted functions and it seems to be 
useful for analysis of soliton type continuous and discrete processes in the N 1 space-time. A close relation 
with the natural generalization of bi- and tri-linear operators into multiple linear operators concludes the paper. 

PACS: 05.45.Yv, 02.30.Jr, 02.30.Gp. 

1 Introduction 

The main goal of this note is the presentation of the role of the addition property (AP), its 
relation to the famous bilinear operator formalism and its universality, since using AP, the quasi- 
periodic processes and soliton ones can be considered in an identical manner. An important 
step is however the generalization of the leading to the bilinear operator standard version of the 
AP to the version which can be linked with multilinear operators. It seems that this generalized 
version can be useful in a case of multidimensional soliton type problems. 

There is opinion that a huge success of the bilinear formalism in the soliton theory can be 
related to the AP for t— functions which appear in the majority of soliton equations 

r (z + w) T (z - w) = ^ We (w) Z, (z) , (1) 

6 

where z = kx + cut e , e e Zl; t : ^ C and We (vf) , Z^ (z) : x Z^ ^ C. 

The essential point here is the factorization of the right hand side of (]l|) , in which functions 
We and Ze depend on w and on z, respectively and exclusively. There are a few version of AP, 
according to (|l]) scheme, [D, and the factorization appears in each one. In apphcations 

to the soliton type equations, the argument z usually depends on space and time, while w plays 
a role of a fixed constant parameter. In a few papers Q it was shown that (|l|) allows in a 
straightforward manner to determine derivatives of logarithms of r— function, which are useful 
in differential version of soliton type equations. For the discrete soliton type equations the form 
(|l|) has a direct and immediate application. 

As it was shown in the cited references, a class containing exponential functions as well as 
the Riemann thcta functions has just the AP according to (|^). 

In order to illustrate an application of AP, we present below two examples: the discrete 
Hirota equation and doubly discrete sine-Gordon equation (dd-sGe). In the hmit, when step 
tends to zero, the first - Hirota equation has a trivial limit, while the second one (dd-sGe) 
becomes a standard sGe. 



2 The Hirota equation 

As an elementary example of the addition property we present the system of dispersion equations 
for the functional (discrete) equations 

or (x -I- h, y, t) t {x — h, y, t) -f 6r (x, y + h,t) t {x,y — h, t) + ct (x, y,t + h)T (x, y,t ~ h) = 
= Cr2(x,y,t), (2) 

where h is the step and a, 6, c, C— are constant. When C = this equation is known as the 
Hirota equation and then its soliton solutions one can be found in pi. Quasiperiodic solutions 



are reported in and another class of solutions in ||^ . In the language of bilinear operator (H) 
can be written as 

[aexp {hD^) + bexp [hDy) + cexp {hDt)] {tot) ^ Ct"^ , (3) 
where bilinear operator (in scalar version) is defined as 

{DxY'' (r o r) := {d^ - d^^)"^ t {x,y,t) t {xi,y,t) \x^^x = (a^)" [t (a; + s,y,t)T{x - s,y,t)] \s=q. 

(4) 



Assuming that r— function argument is z = hx + ly + wt G , equation (|^) can be rewritten 



as 



or (z + kh) r (z - kh) + for (z + Ih) t (z - Ih) + ct (z + w/i) r (z - w/i) = Cr^ (z) , (5) 
which is just suitable for the application of the addition property. We obtain functional equation 

[aW^ (kh) + bW, (Ih) + cW, (wh) - CWe (0)] (z) = 0, (6) 

eezi 

which in case of independent functions (z), £ G Z| leads to the system of algebraic equations 

aW, (kh) + bWe (\h) + cWe (w/i) - CW^ (0) = 0; for each £ £ Z|. (7) 

For the fixed type of solution (which determines the class of functions ) and for a fixed 
step h, equations (0) determine the relations between a, b, c, C and k, 1, w. Thus the dispersion 
equations Q are valid both for soliton and for quasiperiodic processes, and even for processes 
in the form of solitons on the periodic background. 



3 Discrete sine - Gordon equation. 

We consider the functional (discrete-discrete) version of sine - Gordon equation (dd-sGe) in the 
form: 



1 / + U — — U 

7?"n 4 



2 ^^'^ ( 'a ) = sin ( ^ ) , (8) 



where :— u (^ ± h,T ± h) . It is obvious that in the limit h —>■ equation (^) becomes the 
traditional sGe in light cone coordinates 

w,^T = sinu. (9) 

We look for the quasi-periodic solutions (|^) in the form, which is identical (up to constant 
parameters) with solutions of (^) 

e(z+mB) 

"(^'-) = 2*l^-^(^' 

where 9 (z\B) denotes the Riemann theta function of argument z = k^+urG and parametrized 
by the Riemannian matrix Be C^^^; e = [1, IJgZ^, see e.g. Moreover, in order to 

obtain the real solutions we require 9 (z+||i?) = [9 (z\B)]* . 
Since 

9 (z+(k±u)h+mB) 6i(z-|-w±+f IS) , ^ 

u(^ + h,T±h) = 2i In \, ' =2i\n^, ^ \ ' , 11 

' ' 9{z+(k±u)h\B) 9{z + w±\B) ' ^ ' 



and similar relations hold for u — h,T ± h) , after simple manipulations we can rewrite (^) 
in the form 



6* (z + w_) 6* (z - w_) - /i26' (z + w_ + e/2) 6* (z - w_ + e/2) 

6'(z + w+)0(z- w+) ^ ^^'^^ 

6* (z + w_ + e/2) 61 (z - w_ + e/2) - h^9 (z + w_) 6i (z - w_) 
6* (z + w+ + e/2) (z - w+ + e/2) ' 

where w+ = (k + u) h, w_ = (k — u) h. Let us assume that they are constant (with respect 
to ^ and r). Both lead to the same result 

e{z + w_) 9{z- w_) - h^9 (z + w_ + e/2) 6* (z - w_ + e/2) = (C + 1) 6* (z + w+) 6* (z - w+) . 

(13) 

Now the addition property can be applied. The Riemann theta functions do have the 
addition property : 

9 (z + w\B) 9 (z - w\B) ^ ^ We (w) (^z+| |b) , (14) 

for z, w G C^. We coefficients can be expressed by theta-constants but - and it is important 
- do not depend on ^ and r. Then equation ([T^) can be written as 



5: We(w_)02(z+||s)-/^2^ We(w_)0Mz+^)|i?) = (15) 

= (C + 1) ^ We(w+)02(z+||s), 
or as a simple functional equation 

Y [We(w_)-/l2We-e(w_)-(C+l)We(w+)]02|^Z + ||s) =0. (16) 

Since 9^ (z+||_B) labeled by e G ^| form a set of linearly independent functions, finally we 
arrive at the requirement that for any e G Z| 

We (w^) - h^We-e (w_) - (C + 1) We (w+) = 0. (17) 

Equations (|l^) determine of k and u, (and C) and these are linearly related to the prop- 
agation vectors n and angular frequencies w in laboratory coordinate system (z = kx + cot). 
Therefore these equations represent a system of dispersion equations for the discussed dd-sGe. 
Nontrivial solutions of ( p7| ) determines the solutions of starting equation (||) , but for the higher 
g, since then the system is overdetermined, also some and even all elements of the matrix B. 

As (dd-sGe)^(sGe), with /i ^ 0, ( p7| ) also tends to the dispersion equation for standard 
sGe 



Y hUjWe^^j + - (<5e,e ^ cA.fl) = 0. (18) 

id 

where W^^ij := qJ^.q^ (w) |„,=o- In order to prove this statement one can substitute a new 
constant c — —C/h^ .and use the relations We (z) = We (— z), We (0) — Se.o , where S represents 
the Kronecker symbol. 



4 Tri-linear operator 



In order to extend direct methods in the spirit of bilinear operator formahsm on the a broader 
class of equations, the trilinear operators T and T*were introduced by R. Hirota, (c.f. ||, §) 

(T)" {t OT ot) [dz + jdwi + fdw2) t (z) t (wi) r (^2) \w2=wi=z , (19) 
{T*f{TOTOT) := {dz+fd^, +jd^.,)r{z)T{wi)T{w2)U=n„=z , (20) 

where j — exp (i27r/3) . 

In this language e.g. the 5-th order equation of the Lax hierarchy 

U5x + lOuMSa: + 20UxUxx + SOu^Ux + Uf = 0, (21) 

although lacking a bilinear representation, can be written in trilinear form 

{7T^ + 20T^T*^ + 27TxTt) F o F o F = 0, (22) 

where u = 2 (InF)^^ . We will return to this equation in the last paragraph of this note. 
However, first let us try to generalize the concept of bi- and tri-linear operators. 



5 Multiple addition property 

One can introduce multilinear (J-linear) operator by the relation 



(T)" (r o .... o r) := 



J elements 



J-1 



i=0 



(23) 



zj-i = ...=zi=zo=z 



where j — exp {i2-K / J) . 

We are convinced that the effectivity of the bi- and tri- linear operators formalism depends 
on the relevant addition property of the functions to which this formalism is applied. Therefore, 
the fundamental question is which class of functions has the multiple AP. Instead of the class of 
exponential functions appearing in the solutions of soliton type equations we focus our attention 
on the Riemann theta functions constituting a more general class of functions and expressing 
the quasi-periodic solutions. Obviously, exponential functions can be considered as the limiting 
case of the theta functions. 

Following j^] , 1^ , , Q , we adopt the definition of the Riemann theta function as 



I (z|i?) = exp [iTT (2 (z, n) + (n, Bn))] , (24) 



where z e C^, B G C^^^ is the Riemann matrix, (i.e. symmetric with positively defined 
imaginary part) and (z, n) :— Zjrij. 
If z,u('=) eCf, fc = 0, .., J-1, and 



J-1 



5]u(^) = 0, (25) 



fe=o 



one can prove that 



B] = 



(z + u(")|b) 0(z + u(i)|b) 0(z + u(2)|b) ...0(z + u('^-i) 

= ^ I^,(u(«),...,u(-'-i)) ZAz) , (26) 



where 



Z,(z) 

W^e (u("\...,u(^-l)) 

= exp ^i27r < u^"^ 



exp [iir (2 (z, e) + (e, Be))] 9 ( Jz + Be| JB) 



(27) 
(28) 



> 



/ u(o)-u(i)+Se 
u(o) - u(2) + 

V u(o) - u(^-i) + Be 



2B B 
B 2B 

B B 



2B 



B 
B 

2B 



J 



All 6— functions are of order g, except the appearing in (p8[), which is of order (J — I) g. The 
sum in ( ^ ) is over e G i.e. over dimensional vectors whose components are 0,1,..., J— 1, 
and therefore the sum contains elements. Equation (p6|), written here for theta functions, 
is a natural generalization of AP (|^) . The same form has the generalized AP for exponential 
functions defined by 



E (zlB^j = ^ exp in (2 (z, n) + (n, Bn 



(29) 



which appear in the solutions of standard soliton equations. (Observe that the difference between 
( p^ ) and ( p9| ) is only in the number of elements in the sum). It is convenient to assume that 
diagonal elements of matrix B S C^^^ are real. The constraint (BS) can be eliminated easily 
by introducing new parameters w^^^ w'^'^^^^ instead of u*^°^ , u*^.., u*^'^"^) according to the 
relation 



J 



1_ - juC^+i) (1 - Skj-i) - M'^5k,j-i) , k^l, J - 1, 



where 6k j-i is the standard Kronecker symbol. Inversely 



fc-i 



(m) 



J-1 

E 



] w 



(m) 



fc=l,...,J-l, 



J-1 



(30) 

(31) 
(32) 



rn—1 



In the Tabic 1 below we present relations between u^^^ and w^^^ for J — 2 — b. 



J = 


2 


3 


4 






5 






3 = 


-1 


exp (i27r/3) 


i 






exp (i27r/5) 






u(0) = 




w(i) + w(2) 


w(i) + 


W(2) + 




+ w(^) + 


w(3) + 


W(4) 


u(i) - 








j2w(2) - 


f j3.^(3) 




j3-^(3) _ 


f j4w(4) 


U(2) = 






j3w(l) _ 


h jw(2) - 


f j2.^(3) 


j4w(l) + jw(2) + 


^2^(3) _ 


f j3w(4) 


U(3) = 










+ jw'^) 


J^wd) + /w(2) - 


1- jw*^) - 


f j2w(4) 


U(4) = 














F j4w(3) 


+ jw(4) 



Table 1. 



Note that the choice of w parameters is not unique. The set adopted here, gives a full cor- 
respondence with trilinear operators introduced earlier in soliton theory . Since for arbitrary 

integer, J the sum J2k=o 3^ ~ 0' ^^^"^ ^^^^ requirement ( p5| ) is fulfilled for any set w^'^. 

As it was already mentioned, there exist several versions of the addition theorem for theta 
functions. To our knowledge only one form , j^] leads deals to the product of an arbitrary 
number of shifted theta functions as in (p6|), but the r.h.s. is essentially different and unusable 
for our purposes. 

For the fixed J (^^ can be rewritten as 

J 

exp^ [ln^^(z + u(-'')) -ln0(z)] = [0 (z)]"'' ^ (u(i\ u^-^-i)) (z) . (33) 



Differentiating (33) with respect to different components of vectors u'*"'), (fc = 1, J — I) 
we obtain 



Qp+...+q 



■ exp 



j-i 



k=0 
Qp+...+q 



if'^'j - J In 61 (z) 
W,{n^^\...,n 



(34) 



Ze (Z) 



wfiere j = exp (z27r/ J) and u'"' = X]fc=i Cfianging now the derivatives with respect 
u*^*^^ on l.h.s. into derivatives with respect to z, we can find easily the relationship between 
the derivatives of logarithms of theta functions ( with respect to z) and the derivatives of 
functions (with respect to u''^^). All relations become simpler if u^*"') parameters are chosen to 
be zero. 

As an example, the lowest nontrivial derivatives of (u*-^-*) and We (yS'^^ , u*^^^) , respectively, 
(at zero) up to fifth order are reported for J — 2, and 3 below in Table 2 and Table 3. In the 
Appendix we report the lowest nontrivial derivatives (up to 6-th order) also for J — 6. 



J=2 




1/3 

3 X LapL^s + 2L 
15 X LapL^sLi^fi 


+ 4 (^15 X Laf^jsLct, 






Table 2. 



(W^e)„<,, 



{We) 
{We) 



{We) 



2Lal3 
Lap 

4 ^3 X Laf)LjS_ + "^Lap-^sj 

2 ^3 X LapLjS_ + "^Lap-^sj 
ALapL^S + 2 X LajLps + LapjS 

— 2 ^6 X LapL^Scj ~ LaPfSC 

8 (^15 X LapL^Lcj}j + 4 (^15 X L^p^Lcj^j 

~\~'2Lap~^SQ^ 

4 (^15 X ia^L 

+2 ^5 X L gp-ys L^fi + 10 X Lapj^Ls^^ + Lap^scii 

8 ^3 X LapLjsLc_,}j + 12 ^2 X LapLj^Lsfj^ + 
4 X Lapj^Lsp + 4 X L gp-^ f^LsQ + 4Lq^^5L^^-|- 
2 ^6 X LgpLySCti] + Lap-ySCf^ 

4 (9 X LapLysL(;^) + 6 X LasLp^Ljf^ + 

9 X LasLp^^^j^ + 2 ^6 X LapL^sCfj^ + 9 x LapgL^^t^^ 

+ Lgprygl^^ 

Table 3. 



« 7^ j 
i 7^ J 



The remaining derivatives of the order less than 6 - vanish. Lap '■= dz^z^ lii6{z,\B) \z=q 
etc., and we adopted here the shorthand notation which includes all possible permutations with 
respect identically underlined indices: e.g. 3 x LapL^s := LapL^s + La^Lps + L^sLp^. 

Some introductional applications of the above results can be found in Ref. @. 

6 Addition property versus multilinear operators 

The correspondence between the reported here system of dispersion equations and bilinear 
operators is quite obvious. This affinity can be extended even further. For fixed J, let us 
introduce a hierarchy of operators T'-"-' labeled by n = 0, 1, J — 1 

(r("))"'(or)(^-i) := (r("))'" roro...or = (t("))'" [r (z + Uq) ...r (z + u^_i)] luo=...u._,=o 

J-i 

(35) 

where 

,7-1 

t(o) = ^ , (36) 

m=0 

.7-1 J-n 

T(")=a„(o, + rau(™+"-^) + E^""^u(™+"-i), 0^n<J, (37) 

m—J—n+l m—1 

and differentiation relates of course to the indicated components of u*^™^ vectors i.e. Uq™-*. For 
J = 2, 3, 4 we have 





J = 2 




J = 3 




J = 4 








d^m 4 


9u(i) 


9^(0) 4 


9u(l) + 9u(2) 


9u(o) 4 




- 9u(2) + 9u(3) 


r(i) = 


d^m 4 




du(o) 4 


ii9u(i) +i^5u(2) 


du(o) 4 


- ic'ud) - 




r(2) = 










5u(0) 4 




+ j9u(2) + J^9u(3) 


t(3) = 










C^uCO) 4 







It is seen that for J = 2 and J = 3 we have the standard two- and tri-linear operators, 
respectively. However, for J > 3, the operator T^^^ ^ (T*^^^) i.e. T^^^ is not a complex 
conjugate to T^^'. For this reason the operators T^") (n = 1, J — 1) for fixed J, will called 
associated operators. Operator T^^^ is introduced here only for completeness. 

Now, if T— function from ( ^ ) has AP, the question arises how it reflects on the 
We (w(i),...,w("'-i)) functions? 

Using (|T|) and (||) we have d^is) /d^M = ji-^'+P (l + (j-^-i - l) (5o<p<s) and therefore 



= 5^(0) + E ^'"^ui-^-^' + E J'^^uL-^"-) = • (38) 

m—J—n+l m—1 

This means that if r— function has the AP, multilinear operators according to (|35|)-(l37|) 
reduce to the simple differentiation of the (w'^^\ w*^'^^^^) functions with respect to their 
arguments. In the simplest cases of bi- and tri-linear operators this assertion allows immediately 
to write the system of dispersion equation on the basis of bi- or tri-linear approximation. 

As an example, let us note the bi-linear form of Korteweg de Vries equation and tri-linear 
form of the reduction of self-dual Yang Mills equation 



{T^T; + %TlTzT* + grjTt) r o r o r = O, 



(39) 
(40) 



coincide with the relevant dispersion equation systems 



KiUJjWe,WiW, + ^ K.iKjKkKlWe.w,w,WkWi = CWe , (41) 
ij ijkl 

{KiKjKkKlXm - SKiKjKkXlKm) We,WiWjWkWiV^ + ^ X! '^il^j^kWe,w,WjWk = CWe , (42) 
ijklm ijk 

where we assumed that the arguments of r— functions depend Unearly on space and time: 
Zi = KiX + LUit and zi = KiX + Xiy + tOit, respectively. Moreover, in the second equation 
We = We (w^-^), w^^)) L(i)=„,(2)^o, depends on two vectors, designed for typographic reasons 
as w and v. In both cases, C appears as the integration constant and for soliton solutions it 
vanishes, while for quasi-periodic solution it has to be determined as an additional parameter. 
Finally equations ( pi] ) and (^) should hold for any e G Zj, i.e. the first one for e G Z2, and 
the second - for e G 

In conclusion, we expect that the reported here generalized addition property can be useful 
for the analysis of multidimensional soliton equations. 

The Authors acknowledge helpful discussions with Dr. M. Jaworski, Dr. R. Pawlikowski 
and especially are grateful for remarks of Prof. S. Lewandowski. 
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7 Appendix 

The lowest nontrivial derivatives (up to 6-th order) of of (w^^^ w^^^) for J = 6. 









(H^e)„(«„« 






(^e),,(0,,(0,,(j) 

^ ' li^ Uq 






(We)„(>)„(j),,('=) 




« 7^ i 


{Ws) H) H) H) (i) 
^ Ucx ""ff ""y "5 


4 


(3 X LajsLyS + 2Lai3jS^ 






2 


(^3 X La()Lys + 2_L„a,, 


-') 






'^LafjL-fS + 2 X LayLps 


+ Loipjs 








2La/3L~fS + 2 X LayLps 






i, J, fc — di f f event 




3 X La/3 


LyS + Lap-fS 






k, I — different 




-2 (^6 X LapLySi^^ - La 








(W^£)„(*)„(i)„(*)„0)„('=) 




3 X LapLySC^ — 3 X La5Ll3jC 
3 X La^Lp^g — Lap-fSC. 


j ^i,k^i 




-2 X L/SjLaSC - 4 X LacL/SSy - Lafj^sC, 


j j^i,k^i 




+3 X 1 


La-iSLpQ + Lp^gLaC — Lctp^Lsc^ 
2LapLySC ~ Lap-^sc 


different 




-10 X LajS-fLs^ - Lai3jS( 


i,j,k,l,m- 
different 


"" " ' "■' "■- "" 




ir 




+ 4 (ir, X I,,,,-,, 






4 

+ 


^15 X LapLjsL,;^ 

LafJ'fSCfi. 


j + 2 ^5 X La0jsL(;n + 10 X Lap^^Ls^^ 




8 

+ 


^3 X LapLysLc^p^ 
4 X L,,By^,Ls^ + 4 


+ 12 (2 X L^pLy^Lsf,^^ 
LafjjsL^fi + 2 ( 6 X 


+ 4 X Lap-^c^L^ 

L^iSQn 1 + Lap^sCii 




4 

+ 


^9 X LapLysL^fj,^ 
2 ^3 X LapLySCfi " 


+ 6 X LasLp^Ly 

I- 3 X Lap^sLQ^j^ \ 


^ + 9 X LasLpjQfi 

+ 9 X LapsL-y^n + Lap-fSCn 




4 
2 


(3 X LapLjsL^fj,^ 

f 6 X LapLysCfi] 


+ 2 ^12 X LapL-^i^Lsfi 

- 2Lap^jSLi^ii + 6 X La/3 


+ 8 X La<;LpjSn+ 




6 X LaS_Lp^Ly.f^ + 4 ^3 X LapLgi^L^^^ 
+3 X LafiscL-ifj. + 6 X LapspL^Q + 2j 
+2 (3 X LapLrySCfj.) + 6 X LapsL-fQ^j. 


j + 2 (^6 X LapLysL^f,^ 

'^■ySLaPCfi + 2 X LapjsL(^^ 
+ 3 X LaPfj,LyS(^ + Lap-ySC/J. 




6 X LgsLpcL-fu + 2 
+9 X Laps^Lyfj^ + 2 


^9 X LapL~^,sL(^fj_ 

(3 X LapLySCnj 


j + 3 X Lap-fSLQ^ 

+ 8 X LapsLjt^fj. + LaPfSCfj. 




^LapL^gLc^fj, + 2 X i^La-^LpsL^p^ + L 
8 X La^Lp^Lsp + 6 X LapjLs^f^ + 2 

12 X LajS-fC^Sp. + 4 X La-yC,LpSp + -t-o 


CaLi_tpL~i,s + Ly(^Lsfj.Lap^ 
^3 X Lap-fSLi^y^j + 




4:LapLjsL,^ij, + 2 ^4 X LapLji^Lgp^ + 10 x LacLp^L^^ 

+2 {2 X LapLjSCij^ + 13 X Laffj(^Lsp, + 6 x Lap^L^^f^+ 
4 X Laj(^LpSji + Lap-ySCn 




2 ^3 X LapL~fsL(^^ + 12 X La^yLpsL^^ + 6 x Lap^sL^^ 

+ '2-LapLySCti + 8 X La-yL(3&Q^i + 

4 X LafSiLii^^ + 4 X La-isLpQ^ + La/37(5C/i 



